Non-Abelian gauge fields in the gradient expansion: 
generalized Boltzmann and Eilenberger equations 
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We present a microscopic derivation of the generalized Boltzmann and Eilenberger equations in the 
presence of non-Abelian gauges, for the case of a non-relativistic disordered Fermi gas. A unified 
and symmetric treatment of the charge [t/(l)] and spin [S77(2)] degrees of freedom is achieved. 
Within this framework, just as the U(l) Lorentz force generates the Hall effect, so does its SU(2) 
counterpart give rise to the spin Hall effect. Considering elastic and spin-independent disorder we 
obtain diffusion equations for charge and spin densities and show how the interplay between an 
in-plane magnetic field and a time dependent Rashba term generates in-plane charge currents. 
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I. INTRODUCTION 



Spin-charge coupled dynamics in two-dimensional elec- 
tron (hole) gases has been the focus of much theoreti- 
cal and experimental work over the last two decades^. 
Its rich physics belongs to the field of spintronics and 
shows much potential for applications. Thanks to spin- 
orbit coupling all-electrical control of the spin degrees of 
freedom of carriers, as well as magnetic control of the 
charge one, is in principle possible^. Particularly inter- 
esting from this point of view are phenomena like the 
spin Hall effect and the anomalous Hall effect. For a re- 
view of both, see Ref. [H and Refs. 0and|f|, respectively. 
In general terms the theoretical problem at hand is that 
of describing spin-charge coupled transport in a disor- 
dered system. In the semiclassical regime, defined by the 
condition Xp <C I, a Boltzmann-like treatment is sensi- 
ble and expected to provide physical transparency. Here 
Xp is the Fermi wavelength and I a typical length scale 
characterizing the system - say, the mean free path or 
that defining spatial inhomogeneities due to an applied 
field. The Boltzmann equation is a versatile and power- 
ful tool for the description of transport phenomena^, and 
various generalizations to the case in which spin-orbit 
coupling appears have been proposed^—. More general 
Boltzmann-like equations have also been obtained^— . 
In both cases though, much of the physical transparency 
is lost due to a complicated structure of the velocity op- 
erator and of the collision integral. A semiclassical ap- 
proach based on wave packet equation s 14 ' 15 can partially 
circumvent these complications, though it is limited to 
the regime A so t/H 3> 1, with A so the spin-orbit energy 
and t the quasiparticle lifetime. On the other hand it 
was pointed out in different worksi^— that Hamiltoni- 
ans with a linear-in-momentum spin-orbit coupling term 
can be treated in a unified way by introducing SU(2) 
gauge potentials in the model. Taking as an example the 



Rashba Hamiltonian 



Hr 



\-a(p v a x -p x a v ). 

2m 



(1) 



where a is the spin-orbit coupling constant, one can iden- 
tically transform it to 



Hr 



[y + 1 A R a a /2] 
2m 



const. 



(2) 



Here summation over a — X,y,z is implied, 7 is the 
SU(2) coupling constant, and the components of the 
SU(2) vector potential are y(A R )% = -j{Ar) v x = 
2ma, {Ar)% — (Ar)^ — 0. From this point of view a 
different Hamiltonian, say the Dresselhaus one, simply 
corresponds to a different choice of the vector potential. 
An additional advantage of this approach is that it en- 
sures the proper definition of physical quantities like spin 
currents and polarization a 21 ' 23 . More generally, the use 
of the non-Abelian language shows flexibility and poten- 
tial, and has already proven useful in different contexts. 
For example, in Ref. [l9j it was used to predict the exis- 
tence of a "persistent spin helix" in systems with equal 
strength Rashba and Dresselhaus couplings. Such a helix 
was later observed^ and soon after exploited^ 5 -. The au- 
thors of Ref. [13] on the other hand employed it in their 
proposal of a perfect spin-filter based on mesoscopic in- 
terference circuits. Finally, since non-Abelian potentials 
can also be created optically 2 ^—, the range of applica- 
tions of the approach goes beyond systems described by 
Hamiltonians like Indeed, even in solid state systems 
higher-dimensional models such as the one considered in 
Ref. [3(| would fit into the picture^. 

Our goal in the present paper is therefore to put the 
non-Abelian approach in the semiclassical regime on firm 
ground, in order to obtain kinetic equations with a clear 
physical structure and as broad a field of application as 
possible. More precisely, we derive an SU(2) x U(l) 
covariant Boltzmann equation in the framework of the 
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Keldysh^ microscopic formalism. In the covariant ap- 
proach a completely symmetric treatment of the charge 
and spin degrees of freedom is achieved^. Also, we dis- 
cuss the more general Eilenberger equation^— derived 
with the help of the so called ^-integrated Green's func- 
tion technique. The latter allows one to justify the Boltz- 
mann equation in the case when the momentum is not a 
good quantum number due to impurity or other scatter- 
ing, and the notion of particles with a given momentum 
is ill defined. The results obtained hold in the metallic 
regime 6f S/t, with the Fermi energy and H/t the 
level broadening due to disorder, and as long the spin 
splitting due to the [SU(2)] gauge fields is small com- 
pared to the Fermi energy, A so «C e_F, but for arbitrary 
values of A so t/H. We emphasize that, within this ap- 
proach, not only the applied electric and magnetic fields, 
but also the internal spin-orbit induced ones can be po- 
sition and time dependent. 

Our guideline for the present work is the familiar E/(l) 
gauge invariant Boltzmann equation. This reads^ 



Sr + — VR + F-Vp. ]/(T,R,p*) =/[/], (3) 
m ■ 



where the electron distribution function / at time T and 
position R is a function of the gauge invariant kine- 
matic momentum p* = p + eA(T, R) (rather than the 
canonical momentum p), and the Lorentz force F = 
— e[E + (p*/m) A B] appears. The right hand side of 
Eq. contains the collision integral. 

The paper is organized as follows. In Section |TT] we 
start by recalling the quantum derivation of the Boltz- 
mann equation, which allows us to introduce the gen- 
eral formalism in a purposeful way. In Sections Mil the 
generalized Boltzmann and Eilenberger equations are ob- 
tained. In Section IIVI the diffusive regime is discussed 
and spin-charge coupled diffusion equations are derived. 
Finally, Section [Vl shows two example calculations. The 
first involves a study of the Bloch equations in the static 
limit in the presence of in-plane electric and magnetic 
fields, whereas the second is concerned with a novel ef- 
fect, in which an in-plane charge current is generated by 
the interplay of an in-plane magnetic field and a time 
dependent Rashba term. 

We use a system of units where the Planck constant 
h—1 and e = lei. 



II. THE GRADIENT EXPANSION 

The original quantum-mechanical derivation of the 
classical Boltzmann equation by Keldysh 3 - has been ex- 
ploited and extended by many authors, in particular, by 
Langreth^I and Altshuler— . Since it is very instructive, 
we outline the procedure following Ref. dH and consider 
for simplicity's sake the case of free electrons in a per- 
fect lattice. Our aim will be to generalize it to the non- 
Abelian case and to later introduce disorder. The main 



character is the Green function in Keldysh space G 

G A 



(4) 



qR-A are ^fie standard retarded and advanced Green 
functions, whereas G K is the Keldysh Green function 
which carries the statistical information about the oc- 
cupation of the energy spectrum. One starts from the 
left-right subtracted Dyson (quantum kinetic) equation 



-i[G^(l,l')fG(l',2)] =0, 



(5) 



where 1, 1', 2 are generalized coordinates containing space 
and time coordinates as well as spin and Keldysh space 
(and possibly additional) indices. The square brack- 
ets denote the commutator, the symbol "(8)" indicates 
convolution/matrix multiplication over the internal vari- 
ables/indices and 



G,, 



l (l,l') = (id tl 



-tV, 



eA(l)r 



2m 



-e$(l)) 5(1-1') 
(6) 



describes free electrons coupled to an external electro- 
magnetic field. In order to introduce the gradient expan- 
sion we write the Green function in the mixed represen- 
tation in terms of Wigner coordinates 



G(X,p) = / dxe- ipx G(X,x), 



(7) 



where X — ([ti + i2]/2, [xi + X2]/2) is the center of mass 
coordinate and x = (t\ — t 2 , X-i — x 2 ) the relative one, 



X = (T, R) 



(t,r), p = (e,p), px 



-et + p- r. 



Notice that in the presence of both translational symme- 
try with respect to time and space, the dependence on 
X drops out and convolution products as those in Eq.© 
would reduce to simple products in Fourier space. In the 
presence of external fields or in non-equilibrium condi- 
tions, Fourier transforms, as defined in (J7J, of convolu- 
tion products can be systematically expanded in powers 
of derivatives with respect to the center of mass coordi- 
nates. To leading order, the gradient expansion applied 
to Eq.© yields 

- i [Gq 1 f G] « d.G^drG - drG^d.G + 

-V P G 1 • V R G + VrGq 1 • V P G(8) 

Notice that such an expansion is in our case justified by 
the assumption that pf(cf) is the biggest momentum 
(energy) scale of the problem. On the r.h.s. in the above 
both G and Gq 1 are functions of (X,p), with 



Go\X,p) 



[p + eA(X)Y 
2m 



e$(A). 



(9) 



Integrating the Keldysh component of Eq. ((SJ over the 
energy e leads to the l.h.s. of the Boltzmann equation for 
the distribution function 



(10) 
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Note that the above defined quantity is not gauge invari- 
ant. For this reason if one is to obtain a result in the form 
of Eq. ([3]) a shift of the whole equation must formally be 
performed - i.e. one must send p — > p* in Eq. ©. This 
is done in Ref. |37| . though such a shift could also have 
been performed before the gradient expansion. The latter 
is the way followed in Ref. [38j . where the mixed repre- 
sentation is right from the beginning defined in terms of 
the kinematic momentum 

(id 

Unfortunately the simple and convenient concept of a 
"shift" does not work when non-Abelian gauges are con- 
sidered. This is because the nature of the transformation 
([lip is actually geometric, a fact that manifests itself only 
when dealing with non-commuting fields. At its core lies 
the Wilson line Ur(2, 1), which is the exponential of the 
line integral of the gauge potential along the curve T go- 
ing from 1 to 2 (see [391 ]) 

C/ r (2,l) = Ve- iri * dyA{v) . (12) 

Here the symbol V stands for path-ordering along P, 
whereas 77 is a general coupling constant - in the U(l) 
case it reduces to e. In the spirit of the gradient expan- 
sion the integral (fT2")l is evaluated for small values of the 
relative coordinate x, and it is thus reasonable to pick T 
as the straight line connecting the two points. For the 
[7(1) gauge it is seen that 

U r (2, 1) w e - ieA{x)x (13) 

which is precisely the phase factor appearing in Eq. (1111) . 
In other words, the "shifting" of Eq. (JS| should properly 
be seen as the transformation 

[0^(1, 1')? 6(1', 2)] -> 
U r (X,l) [Go l (l,l') ? G(l', 2)] U r -(2,X), (14) 

where T(r') is a straight line from 1(A) to X(2)22 and 
matrix multiplication over the internal indices between 
Ur, Ur> and the commutator is implied. 

With this hindsight about the nature of the "shift" 
required to obtain Eq. @, it is possible to generalize 
the construction to the non-Abelian case. A general 
gauge transformation is defined as a local rotation of the 
second-quantized annihilation fcrmionic field "0 

^(1) = V(1M1), V(l)Vt(l) = l. (15) 

For the gauge potential one has 

V A'(1) = V(l) foA(l) + idx] V\l). (16) 

Notice that this is now a tensor with both real space and 
gauge indices 

77,4(1) = (e* + 7*^72, eA + 7^t a t a /2), (17) 



where we found convenient to separate the Abelian cou- 
pling constant e from its non-Abelian counterpart 7. The 
non-Abelian scalar potential 'yvE ,a i°/2 describes a Zeeman 
term - i.e. of the kind b • s. Here s is the spin of the 
carriers, whereas b could be an applied magnetic field 
or, in a ferromagnet, the exchange field due its mag- 
netization. The t a /2's are the generators of the given 
symmetry group, which in the SU (2) case become the 
Pauli matrices, t a = a a ,a = x,y,z. In the following 
boldfaced quantities will indicate vectors in real space, 
whereas the presence of italics will denote a gauge struc- 
ture - which, as in the above, will sometimes be written 
down explicitly. A sum over repeated indices is always 
implied. Since a Wilson line transforms covariantly, i.e. 
£/f(2,l) = y(2)[/r(2, 1)0(1), it is possible to define a 
Green's function 

G{l,2) = U r (X,l)G(l,2)U r (2,X) (18) 

which is locally covariant, i.e. 

G'(l, 2) = V(X)G(1, 2)V t (X). (19) 
In terms of G K we can define a distribution function 

/(^P)4(l + /^(A,,)) (20) 

which will be the natural generalization of f{X, p*) from 
Eq. ([3]). The procedure is then clear: 

1. transform the kinetic equation according to 
Eq. HUD; 

2. expand the Wilson lines (see below); 

3. perform a gradient expansion and write everything 
in terms of G K (X,p); 

4. integrate over the energy e to obtain the Boltzmann 
equation, or over £ = p 2 /2m — [i to end up with the 
Eilenberger equation. 

Postponing the discussion of the last point to the next 
Section, we now consider the general expression 

Go 1 (X,p) = e-H(X,p) 

[ P + e A(A)+7^ a (A)tV2] 2 
2m 

+e<P{X)+jf a (X)t a /2. (21) 
In the Rashba model, Eq. ©, one for example identifies 

The procedure outlined above (points 1.-3.) leads to 

a locally covariant equation for G accurate to order 
[(dxdp)(Ad p ), (Ad p ) 2 } (see^): in the mixed representa- 
tion language we have formally two expansion parame- 
ters, dxdp <C 1 - the standard gradient expansion one 
- and Ad p < 1 - coming from the gauge fields. In the 
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SU(2) case the latter corresponds to the physical assump- 
tion that the spin-orbit energy be small compared to the 
Fermi one, A so /ep <C 1. Even though our treatment is 
valid for any non-Abelian gauge, we now pick the SU{2) 
gauge for definiteness' sake. In this case steps 1.-3. lead 
to 

(Bt + ~ ■ Vr - ■ {[eE + 7 £] d e , . } + 



1 



{^•V p , .})<? = 0,(23) 



where the symbol { . , . } denotes the anticommutator. 
The covariant ("wavy") derivatives are 



d T = d T - »7 [*,.], V R = V R + vy [A, . ] , 
whereas the generalized Lorentz force reads 



E- 



p AB 



1 



(24) 



(25) 



U(l) 



SU{2) 



The fields are given as usual in terms of the field tensor 
F, but this has now an SU{2) x f/(l) structure 



Ei = = -d T A 



V fl4 $ (26) 

£t = K = -9tA? - V R , * a + z 7 [¥, Ai] a (27) 

Bi = le ijk ¥%, Bf = ^e ijk W% (28) 

¥ a ]k = V Rj A k -V Rk A 3 (29) 

F« fc = V Rj Al - V Rk A« + i 1 [A„A k } a . (30) 



Note that in order to obtain Eq. (|23[) it is sufficient to 
expand the Wilson lines to first oder in x, e.g. 



U T {X,1) wl + t/il 



IX 



(31) 



This is not true for a general convolution of the kind 
[F(l, V) f G(l', 3)], with F(l, 1') a function with a more 
complicated structure than that of Gq (1,1'). Such a 
case would require a second order expansion, e.g. 

Ur(X, 1) « 1 + ^| + - rj 2 A 2 ^-, (32) 

and would lead to a rather more complicated equation. 

To complete our preparatory work for the derivation of 
the kinetic Boltzmann or Eilenberger equations, we need 
to introduce the effect of disorder. Within the Keldysh 
formalism this is done by the addition of a self-energy 
contribution on the r.h.s. of Eq. <j5j> 



[£(1,2)? 6(2, 3)] , 



(33) 



which can be manipulated just as the "free" (Gq 1 ) term. 
In spin-orbit coupled systems the presence of disorder 
can have a number of interesting effects. Indeed, phe- 
nomena like the spin Hall effect, anomalous Hall effect or 
related ones can have both an intrinsic and an extrinsic 



origin . This depends on whether they arise from fields 
due to the band or device structure, or from those gen- 
erated by impurities. In the latter case skew-scattering 
and side-jump contributions to the dynamics appea r 42 ' 43 . 
For a discussion of these issues see^i^— . In the following 
we limit ourselves to the treatment of intrinsic effects in 
the presence of spin independent disorder. We consider 
elastic scattering with probability W — W(p — p') and 
quasiparticle lifetime r _1 = 27rVJ p , 5(e p — e p i)W(p — p'). 
In the Born approximation, the disorder self-energy in 
the mixed representation reads 

S(X,p,e)=^W(p-p')G(X,p',e). (34) 
p' 

From Eq. (|34l) one obtains that the locally covariant self- 
energy S is 

E(X,p,e)=^W(p-p')G(X,p',e), (35) 
p' 

which in turn implies 



I[G] = -i 



E, G 



(36) 



Note that for a leading order description of the cou- 
pling between spin [SU(2)} and charge [f (1)], corrections 
O(Adp) in the collision integral are enough. Notice also 
that, whereas G is peaked at the different folds of the 
spin-split Fermi surface, the peaks of G arc "shifted" and 
thus located on the Fermi surface in the absence of spin- 
orbit coupling. 



III. THE BOLTZMANN AND EILENBERGER 
EQUATIONS 

The question of whether to integrate the locally covari- 
ant kinetic equation with respect to e or to £ = p 2 /2m — fi 
depends on the physical situation. If the spectral density 
i(G R — G A ) is not "(5-like" as a function of e the energy 
integration is formally impracticable. The ^-integration 
on the other hand is capable of justifying a Boltzmann- 
like approach even when the first approach fails or looks 
severely limite d 34 ' 35 . For the case considered of a degen- 
erate gas of free electrons colliding elastically with impu- 
rities both procedures are viable, provided the condition 
tpT 3> 1 holds - since, as mentioned before, all quantities 
appearing in Eq. Q36p are peaked at the A so = Fermi 
surface. 



A. Boltzmann (e-integration) 

Energy integration of the Keldysh component of 
Eqs. (f2"3"| and (|36l) yields a Boltzmann-like kinetic equa- 
tion for the 2x2 matrix distribution function f(X, p) 

d T + £ • V R + i {T ■ V p , . }] f(X, p) = /[/] (37) 
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with the covariant derivatives and the generalized 
Lorentz force T defined respectively as in Eqs. (j2~4")l and 
(|25| . and where the collision integral reads 



I[f] = -2nJ2W(p-p')6(e p -e p >) [f(X,p) - f(X,p')] . 
p' 



(38) 

Notice that Eqs. (|37 ]) -(p8 |) are formally valid both in two 
and three dimensions. However, since the physical sys- 
tem we have in mind is a two-dimensional electron gas, 
from now on we restrict ourselves to two dimensions. 
Observable properties are conveniently expressed via the 
matrix density, p, and current, JT, 



d p 

]2 



/(*,P) 



(39) 



*( X ) = / 7TT2 Z /(^P)' ( 4 °) 
J (2tt) 2 m 

which obey the following generalized continuity equation, 

d T p(X) + V R • J{X) = 0, (41) 

derived by integrating Eq. (|37p over the momentum. Ob- 
servables like the the particle and spin densities, n and 
s a , and the particle and spin currents, j° and j a , can be 
evaluated as 



n(X) = Tr [p{X)] 
S a (X) = \Tr[a a p(X)} 
j°(X)=Tr[J] 
■f(X)= l -^[a a J]. 



(42) 
(43) 
(44) 
(45) 



One can check that these expressions agree with their 
microscopic definition o 18 ' 21 . 

Eq. (f3"7| is the first main result of the paper. Though 
the idea of rewriting spin-orbit interaction in terms of 
non-Abelian gauge fields is no novelty, we are not aware of 
a Boltzmann formulation in the above form. Whereas in 
Refs. [8jj9[ the collision integral and the velocity are non- 
diagonal in the charge-spin indices, here their structure 
is simpler. The gauge fields appear only in the covariant 
derivatives, describing precession of the spins around the 
external magnetic field and the internal spin-orbit one, 
and in the generalized Lorentz force, which couples the 
spin and charge channels. 



B. Eilenberger (^-integration) 

The integration over £ of Eqs. (|2"3"|) and (f3"fj)) yields in 
two dimensions the Eilenberger equation 



1 



(&t + v F p • V R - ^d e 



P(e) 



(eE + 7 £) 



1 



+^{F(p F ,v)-[-P + <pd v ,], ■})g K 

= -2^o / - J) [g K (<p) ~ g K ^')](m 



where p = (cos (p, sin ip), (p — (— sine/?, cos ip), W(ip — ip') 
is the scattering amplitude at the Fermi surface, and g K 
is the Keldysh component of the covariant quasiclassical 
Green's function 



dCG(X,ip,e,0- 



(47) 



Notice that the energy derivative d t acts on the whole 
anticommutator, i.e. on g K too. Just as in the Boltz- 
mann case, and as opposed to what happens in the 
literatur o 12 ' 13 , the velocity and the collision integral have 
here a simple diagonal structure, whereas the gauge fields 
appear only in the covariant derivatives and force terms. 
The collision integral will be extensively discussed in Ap- 
pendix lA 21 to make an explicit comparison with Ref. |12j | 
possible. Integration of Eq. (|46|) over the energy and the 
angle leads again to the continuity equation (|4T|) . this 
time with densities and currents expressed in terms of g 

P{X) = J M~9 K ) (48) 

J(X) = -^- J dev F (pg K ), (49) 

where (...) denotes the angular average. Recall that when 
expressing physical quantities in terms of the standard 
quasiclassical Green's function g — i/ir J d£G equilib- 
rium high- energy contributions are misse d 34 ! 36 . For in- 
stance, Eq. (|48]) for the particle density when only U(l) 
fields are present would be written, in terms of g K , as 



/>..V)=-^ / A<[^) - V,„<I»(A-) 



(50) 



with the second term due to the scalar potential originat- 
ing from the high-energy part. A virtue of the present 
formulation is that such contributions are by construc- 
tion included in the covariant g. Moreover notice that, 
whereas in the presence of spin-orbit coupling the usual 
normalization condition g 2 — 1 is modified and be- 
comes momentum-dependenlii 3 -, in the covariant formu- 
lation g 2 — 1 holds - see Appendix [A] The normal- 
ization condition is established by direct calculation "at 
infinity" , i.e. where, far from the perturbed region, the 
Green's function reduces to its equilibrium form. It plays 
the role of a boundary condition imposed on Eq. (|46l) , and 
thus defines its solution uniquel y 48 i 49 . In the presence of 
interfaces between different regions wave functions have 
to be matched, and this can be translated into a condi- 
tion to be fulfilled by the quasiclassical Green's function 
on either side of the interface d 5 ^ 51 . Recently some very 
general such boundary conditions for multiband systems 
were obtained 5 - 2 -, though valid only as long as the spin 
and charge channels are decoupled. When this is not 
anymore the case, things are complicated by the momen- 
tum dependence of the normalization and, as far as we 
are aware of, beyond the present treatment of bound- 
aries. The covariant formulation in terms of g suggests 
however the possibility for a non-trivial extension of the 
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known boundary conditions to the case in which spin and 
charge channels are coupled, precisely because of the sim- 
ple normalization of g. 



IV. THE DIFFUSIVE REGIME 

Our goal in this Section is the discussion of spin- 
charge coupled dynamics in the diffusive regime. For- 
mally, the "non-Abelian" Boltzmann equation ([37]) can 
be solved just as in the U(l) case. We expand the 
angular dependence of the distribution / in harmonics, 
/ = (/} + 2p • f + ... and use the expansion in Eq. (|37|) 
to obtain an explicit expression for f , 



with 



2m 

Ttr 



Ttr 

2 



V P , {/)}} + 



-y(p{^-V p ,(2p-f)}) 

= fdiff + fdrift + feall- 

In the above Tt r is the usual transport time 



(51) 



— = 2vrJV f ^W(<p - <ff)[l - cos(^ - <p% (52) 

Ttr J 27T 

which depends on the energy £ = e p —e F through the scat- 
tering probability W (<p - <p' ) = W{£ , <p - <p') |f =f . The 
diffusion term i^is is related to the (covariant) derivative 
of the angular average of the distribution function, i.e. to 
the derivative of the charge and spin densities. The drift 
term fdrift arises from the second term on the r.h.s. of 
Eq. (f5Tj) . in which only the "electric" part of the Lorentz 
force, i.e. — [eE + 7 £], contributes. The Hall component 
fkaii comes instead from the third term on the r.h.s. of 
Eq. (|51[) . and is due to the "magnetic" part of the Lorentz 
force, -p A [eB + 7 B] /m. Using Eq. flST} into Eq. ([10]) 
one finally has 



J 



dp P 

(2tt) 2 m 

d 2 p 



[</} + 2p-f + .. 



/ (2tt) 2 m 

J diffusion + J drift 



2p • [fdiff + fdrift + feall] 
J Hall- 



(53) 



The drift current is straightforwardly computed 

Adrift = ~ Mm), eE + 7 £} , u{p) = -N D{p) (54) 

where No is the density of states, D(/j,) the energy de- 
pendent diffusion constant and /x the (spin dependent) 
electrochemical potential. Since we assume fields that 
are small compared to the Fermi energy, it is often suf- 
ficient to replace D(/i) by its value at the Fermi energy 
and in the absence of the U(l) and SU{2) fields. In the 
examples we discuss it will be important to go one step 
beyond this simple approximation, for which we obtain 



D(jjl) w D(e F ) + d^D (p - N e F )/N Q 



(55) 



D(e F ) 



V F T tr 



0^ = ^(l+7o/2), (56) 
m 



/o mv F a 
7o/2 = -dt-Ttr- 

Ttr 



(57) 



The factor 70 is defined to make direct contact to 
Ref. [53j]. Notice that due to the expansion in (1551 the 
diffusion constant D(p) becomes a spin dependent ob- 
ject, 



D{n) = D° + D a a a 



(58) 



with D° ps D(e F ) and D a ps d i Ds a /N . 

The calculation of the diffusion current is slightly more 
involved: the momentum integration is delicate, since the 
integrand has a nontrivial matrix structure and is out of 
equilibrium. In order to "extract" such a structure we 
first write 



J di 



ffusion 



(59) 



thus defining a diffusion constant T> which is now a ma- 
trix. Extending the Einstein relation to the present non- 
Abelian case will give T> an explicit form. At equilibrium 
one has 



iV [e$ + 7*] + N e F , 



(60) 



and as - again, at equilibrium - the diffusion current 
balances out the drift one 

J drift = —Ja iffusion 

= -^{V,N [eE + 1 S}} (61) 



there follows 



V = D(n) (62) 



fhaii = ^{eB + 7 BAf}, 



from which we find 



as to be expected. 

The Hall term fnaii can be obtained from the equation 
implicit in Eq. (|5T1) 



(63) 

(64) 

(65) 
(66) 



with 



<7 H aii = — {eB + jB ^ { rtr (/i), J}} . 
Am 



Ttr(fJ-) = T tr + d ( T tr (p - N e F )/N 

7o T tr p- N e F 



Ttr 



2 mvp 



iVn 



To be more explicit we give the expressions for the 
particle current j and spin current j a 



V out (0 



j° = -D\\7n + 2eN E) -2D a ([Vs} a - 

_^j0 AB _2Z*Lj» AB a 

m m 

■f = --D a (Vn + 2pJV e) - d([Vs] u 
— ——j a A B — -7-^-j A B a . 



7^0 ca 



(67) 



(68) 



Here we included the Hall current only in the leading 
approximation, i.e. Tj r (/i) ~ T"tr(eF)- The diffusion 
equations for charge and spin are obtained by inserting 
Eqs. (|67| and (|68|) into the continuity equation (|41~j) . 



TWO EXAMPLES 



VinCO 




V out (t) 



A. Effect of an in-plane magnetic field 

As a first simple example that shows how the formalism 
works, we obtain and solve the Bloch equations for a 
Rashba 2DEG driven by an electric field along x and in 
the presence of an in-plane Zeeman field along x. This is 
the same geometry considered in Refs. [UIHS andlH3|. 

The Z7(l) fields read 

E=(£,0,0), B = (69) 

while, since the Zeeman field enters the Hamiltonian 
through the scalar potential j^> x a x /2 = b x <r x /2, the 
SU(2) ones are 

7 £ = 2mab x (<r z /2, 0, 0), jB = -(2ma) 2 (0, 0, a z /2). 

(70) 

From the expressions for the currents derived in the pre- 
vious section, one obtains in the homogeneous limit a set 
of Bloch equations which generalizes those appearing in 
Ref. [13] to the case of angle-dependent scattering - but 
in the absence of extrinsic effects -, namely 

s = -f [s - bN /2 + eaT tr N z A E] + 

-[b-2ear ir (l + 7o/2)zAE] As. (71) 

Here T — 1/tdp diag(l, 1, 2) is the relaxation matrix, 
with 1/tqp = (2ma) 2 D the Dyakonov-Perel relaxation 
rate. Notice that the electric field in the first and in the 
second term on the r.h.s. of Eq. (71) has a different 
origin. While the first term is traced back to the (spin) 
Hall current and therefore to the SU(2) magnetic field, 
the second term can be traced back to the drift current. 
The factor 70 that appears due to the energy dependence 
of the scattering time has an important impact on the 
static solution of the Bloch equations. When 70 = we 
find 



FIG. 1. (color online) The Rashba spin-orbit coupling con- 
stant a is made time dependent by applying a time depen- 
dent gate potential Vi„(t). The light (blue) area represents a 
two-dimensional electron gas inside a heterostructure. When 
an in-plane magnetic field along x, b x , is also switched on, 
a charge current j® flowing along y and proportional to a is 
generated, its actual sign depending on the sign of a. The 
induced voltage drop in the transverse direction V ut(t) can 
then be used to measure the strenght of the Rashba interac- 
tion. 



i.e. the effects of the Zeeman and the electric field on the 
spin polarization are simply additive. This is not any- 
more the case if 70 7^ 0, in which case we find in the 
limit of weak electric and magnetic fields (in our geome- 
try both in x-direction) 



-eaT tr N E, 



-70- 



eaT tr N Eb x 



4{2ma) 2 D ' 
(73) 

that is, in-plane fields generate an out-of-plane spin 
polarizatio n 5355 .In the above s x = b x No/2. 



B. Charge current from time dependent spin-orbit 
coupling 

The Rashba spin-orbit coupling constant a arises from 
the potential confining the 2DEG and is thus tunable by 
a gate voltage: if the latter is time dependent, so is the 
former. Let us then consider the Rashba Hamiltonian for 
a time-dependent Rashba parameter, a —> a(T). In the 
non-Abelian language this means that the SU(2) vector 
potential becomes time dependent, and therefore a spin 
dependent electric field is generated. Explicitly we have 



7 £ = 2ma (0-72,-0-72,0) 



s = bA^o/2 - eaT tr N z A E 



(72) 



7B 



{2maf (0,0,ct72), 



(74) 
(75) 
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with a = dra- The SU(2) electric field leads to the 
appearance of in-plane spin currents, as discussed in 
Ref. (56{. However it does not generate a charge current, 
since it acts with opposite sign on particles with different 
spin: the net field obtained after averaging over all par- 
ticles is zero. This is not anymore the case if a magnetic 
field is also present. Say the latter points in x-direction, 
then a nonzero average electric field in the y direction 
appears, given by j(£ y ) = —ma(<7 x ) — —2mas x /n 
(here (...) denotes the average over all particles). We 
then expect a particle current in y-direction of the order 
j y = — 2DNo , y(£y). We now make the argument quanti- 
tative. Let us apply an in-plane Zeeman field along x, as 
shown in Fig. [TJ Then the SU (2) electric and magnetic 
fields are 

1 £ = 2m (aa y /2 + ab x a z /2, -aa x /2, 0) , (76) 



7 



(2ma) 2 (0,0,a72) 



(77) 



Note that the structure of the Bloch equations (|7T|) is not 
modified, so that the stationary spin density is 



a x — B x B y - 



0, s z = 0, 



(78) 



with s x q = N b x /2 as before. As expected, the SU{2) 
electric field generates a particle current flowing along y, 



No 
2 

T tr aN b x (l+ l0 /2), 



fy = -2D X ^j£* 



(79) 



having used D x — T tr /m(l + 'jo/2)s x . Finally, for a gen- 
eral direction of the in-plane magnetic field b the charge 
current is given by 



j° = T tr N a(l + 7o /2)zAb. 



(80) 



Such an effect could provide an alternative way of es- 
timating the strenght of the Rashba interaction, since 
other spin-orbit mechanisms would not gain any time de- 
pendency from a modulated confining potential. 



VI. CONCLUSIONS 

We showed how to microscopically derive the general- 
ized Boltzmann and Eilcnberger equations in the pres- 
ence of non-Abelian gauge fields. In the SU{2) case 
such equations can be used to describe spin-charge cou- 
pled dynamics in two-dimensional systems whose Hamil- 
tonians include linear-in-momentum spin-orbit coupling 
terms. All degrees of freedom are treated symmetrically 
and the proper identification of the physical quantities 
follows naturally from the form of the continuity equa- 
tion. Considering elastic disorder, we obtained results 
which hold as long as £ > 1/r, A so and for arbitrary val- 
ues of A so t. In particular, we showed that by using the 
covariant quasiclassical Green function, the collision in- 
tegral in the kinetic equation is not affected by the gauge 
fields, which only appear to modify the hydrodynamic 



derivative. We expect that this nice disentanglement of 
gauge fields and disorder effects in the Boltzmann and 
Eilenberger equations may prove very useful when con- 
sidering quantum correction o 57 i 58 . We also expect the 
approach to allow for a generalization of the boundary 
conditions for the Eilenberger equation to the case in 
which spin and charge channels are coupled. When dis- 
cussing the diffusive regime, we first obtained Bloch-like 
equations for the spin and charge, and then exploited 
them to predict a novel effect. Finally, we note that 
by making the non-Abelian coupling constant momen- 
tum dependent, 7 — > 7(p), it may be possible to extend 
the present formalism to include Hamiltonians with more 
general forms of spin-orbit interaction - i.e. not limited 
to being linear in momentum. 

We thank U. Eckern, J. Rammer and M. Dzierzawa for 
discussions. This work was supported by the Deutsche 
Forschungsgemeinschaft through SFB 484 and SPP 1285 
and partially supported by EU through PITN-GA-2009- 
234970. 



Appendix A: Quasiclassical normalization condition 
and the collision integral 

Consider a quantity F(l,2) which is non-locally 
covariant, i.e. which under the gauge transforma- 
tion Eq. (fT5")) transforms according to F'(l,2) = 
V(1).F(1, 2)V <! (2). Its locally covariant counterpart reads 
F(l,2) = U r {X,l)F(l,2)U r (2,X), where X = ([h + 
t2\/2, [xi +x 2 ]/2). In Wigner coordinates, up to 0{Ad p ) 
accuracy, one has 



F = F--{Ad p ,F}. 
We define the ^-integrated functions 



f{e,<p,X) = ± J d£F(p,e,X) 



Ke,<p,X) 



d£F(p,e,X) 



(AI) 



(A2) 
(A3) 



Let us start by assuming for simplicity A = (0, A. a a a /2) 
- that is, we have neither electric nor magnetic fields, 
only spin-orbit coupling - and setting the SU(2) coupling 
constant to one, 7=1. Moreover, the functions F, F are 
assumed to be peaked at the Fermi surface £ = or in 
its vicinity. Thus, by ^-integrating Eq. (|A1[) by parts one 
has 



/ = / + 



1 f Ap A-ip 



Pf 



Pf 



d v ,f 



(A4) 



with pf the Fermi momentum in the absence of spin-orbit 
coupling. The presence of a U(l) vector potential can be 
handled just the same way, whereas the inclusion of the 
scalar potentials e<& + ^ a a a /2 is trivial and amounts to 
a shift of the energy argument of / 



/ = /-^{(e$ + *)9 £ ,/} 



(A5) 
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Therefore, in the presence of a general 4-potential A 
(e$ + eA + A a a a /2), one has 



f = f-\{(e$ + *)d e ,f} + 

1 f (eA + A)-p (eA + A) ■ <p , 



d v ,fUA6) 



2 i p F p F 
We now use Eqs. (jAl|) - (|A6|l to show 

• how the ^-integrated Green's functions g and g are 
related, and what this implies for the latter's nor- 
malization; 

• that the collision integral Eq. (|5r?)) is equivalent to 
the one appearing in [12j. 

1. About q and q 



(A7) 



Take F = G. From Eq. (|A4j) one obtains 



fJ 



A p A ■ (p 



PF 



PF 



where we did not write down explicitly all the dependen- 
cies, since no confusion should arise. Direct calculations 
sho 1 



43 



R,A 



± 1 



A- p 

PF 



4 = 2tanh(e/2T)(l 



A p 

PF 



(A8) 
(A9) 



T being the temperature and the result for the Keldysh 
component being valid at equilibrium. It follows that 
to order |«4.|/p,p - i.e. A so /e_F - the locally covariant 
^-integrated Green's function has no SU(2) (spin) struc- 
ture 



= ±1 



gf 9 = 2tanh(e/2T) 



(A10) 
(All) 



and satisfies the standard normalization condition g 2 = 
1. Recall this has the meaning of a boundary condition 
satisfied by g, and so is not affected by the introduction 
of driving electromagnetic [U(l)] fields^ or of a Zeeman 
term. Indeed, we saw that including the scalar potentials 
e$ and \I/ simply "shifts" the energy argument of g, 



7.K 



2tanh(e/2T). (A12) 



2. The collision integral 

Take F = -«[£, G] K = C, and soF = -i[S, G] K = C. 
The ^-integration delivers 

c=~[(K)g K -(K~g K )], (A13) 



with the kernel K((p — ip') = 2ttNqtW(ip — (p'), and where 
(...) is shorthand for angular average. One then uses 
the inverse of Eq. (|A6|) to calculate the corresponding 
expression in the standard ("non-tilde") language. We 
consider separately the effects of spin-orbit (A) and of a 
Zeeman field (vP), since they add linearly. 
First, spin-orbit. Starting from 



1 

C = C -2 



PF 



PF 



-[g-(Kg)} 



1 

27 



A ■ p A - if 



PF 



PF 



d v ,g-(Kg)}, (A14) 



the translation from g to g is done by means of Eq. (|A4I) . 
The calculation is easy but some care is needed, so this 
is done step by step. First recall that 
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A ■ p A- <p 



PF PF 

so that Eq. (|A14[) becomes 

c = — [g - (K g)] 4 



d Vl g 



9, 



(A15) 



1 

27 



A p A ■ tp 



PF 



PF 



d v , (Kg) 



(A16) 



Then consider the (Kg) term, where K = K((p — ip') 

dtp' 
27 

i r ^ 



(Kg) 



-9(<P') 



+ 



2tt 



2 

(Kg) + 

w<V 

2 J 2n 

A-<p' 



A ■ p' T , A-cp' 



-K 



PF 



PF 



K,d v ,g(ip') 



A-p' t A-(d v .tp') 



PF 



PF 

[d V 'K],g(ip') 



PF 



K+ 



having performed a partial integration and used that 



g{<p = 0) 
reads 



g(ip = 2tt), d^ip' 



T 
1 

27 

1 

27 



(Kg)] 



l 

27 

dip' {A - 0' 

^7 
A- tp 



PF 



PF 

d v ,(Kg) 



p'. This way Eq.(|M6| 

^(Kg) 

PF 

,K>K]g(ip') 



(A17) 



Now work on the last term. Recall the assumption that 
the scattering amplitude depend only on the momentum 
transfer, i.e. K(p, p') = K(p — p'). This implies 



JLq k = -^pd^K - 
PF m 



PF , <f „ 
— P ot H d v > 

m pf 



K. (A18) 
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From the last term of Eq. (|A17|) one therefore has 

Substitution back into Eq. (|Al7j) gives 
c = -~[g-(Kg)} + 

+ 5 L ({ — ,3 ( K(,\). (A20) 

This expression can also be obtained by a direct £- 
integration of the collision integral in the standard ( "non- 
tilde") language. It agrees with the one appearing in 
12] for the case of parabolic bands when one identifies 
b • er/2 = A. ■ p/m, b being the internal spin-orbit field 
in the language of Ref. (l2{ . Besides the first two terms, 
in which no spin-orbit contribution appears, the third 
term corresponds to corrections due to the spin depen- 
dent density of states, whereas the fourth and fifth arise 
from the energy dependence of the scattering amplitude. 
In the notation of (TUl- the former corrections are due 
to M d and the latter to - notice that some of the 
M d and M™ contributions cancel each other because of 
Eq. (TAl8|) . 

Let us now consider a Zeeman field described by the 
scalar potential \&. Shifting from covariant to non- 
covariant quantities is done according to Eq. (| A5|) and 
its inverse. Notice that the kernel K is a function of £ 
evaluated at £ = e. The energy e is actually sent to zero 
during the £— integration, though since we now have to 



shift back to the non-covariant language it is better to 
explicitly keep track of this dependency. At the end it 
will be as usual £ = e — >• 0. From the inverse of Eq. (IA5|) 

c(e)=c(e) + ±{*,d e c(e)}. (A21) 
The first term on the r.h.s. gives 
c=--[(K)g-(Kg)} + 

T 

-±-[(K){9,g'}-(K{9 )g '})], (A22) 
whereas the second one leads to 

d £ c = -- [{K')~g + (K)g' - (K'g) (Kg')] (A23) 

T 

where K' = d e K,g' = d e g. Plugging both expressions 
into Eq. (|A21I) one has 

c(e)=--[(K)g-(Kg)} + 

T 

+-L[(K){*,g'}-(K{*,g>})} + 

--L[(K'){*,g} + (K){*,g>} + 
" -(K'{9,g})-(K {%</})] 
= -\[{K)g-(Kg)] + 

-!-{*,{K>}g-(K'g}} (A24) 

with g = g(e),g' = d e g and K,K' are evaluated at the 
Fermi surface £ = e = 0. The full expression for the 
scattering kernel in the presence of spin-orbit coupling 
and a Zeeman field is given by the sum of Eq. (|A20|) and 
Eq. (|A24|) . It leads to results in agreement with Ref. [53| . 
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